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We consider a dynamic method, based on synchronization and adaptive control, to estimate 
unknown parameters of a nonlinear dynamical system from a given scalar chaotic time series. We 
present an important extension of the method when time series of a scalar function of the variables 
of the underlying dynamical system is given. We find that it is possible to obtain synchronization 
as well as parameter estimation using such a time series. We then consider a general quadratic flow 
in three dimensions and discuss applicability of our method of parameter estimation in this case. 
In practical situations one expects only a finite time series of a system variable to be known. We 
show that the finite time series can be repeatedly used to estimate unknown parameters with an 
accuracy which improves and then saturates to a constant value with repeated use of the time series. 
Finally we suggest an important application of the parameter estimation method. We propose that 
the method can be used to confirm the correctness of a trial function modeling an external unknown 
perturbation to a known system. We show that our method produces exact synchronization with 
the given time series only when the trial function has a form identical to that of the perturbation. 
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I. INTRODUCTION 

An experimental observation often consists of reading 
a time series output from a dynamical system. Such a 
time series can contain information about the number as 
well as the form of the functions governing the evolution 
of the system variables including nonlinearities (if any) 
and the parameters jlj . The estimation of parameter val- 
ues from a given chaotic scalar time series of a nonlinear 
system is the topic of our interest here. 

We have recently given a method to dynamically esti- 
mate unknown parameters from the chaotic time series of 
a single phase space variable when the system equations 
are known [0 . The method is based on a combination of 
synchronization |^-|^] and adaptive control similar to 
that used by John and Amritkar 

The problem of parameter estimation in nonlinear dy- 
namics has been considered earlier. Parlitz, Junge and 
Kocarev have given a static method H based on min- 
imization while Parlitz has developed a method based 
on auto-synchronization [ fl0|| . Unlike our method, auto- 
synchronization method requires an ansatz for the pa- 
rameter control loop and gives slower convergences in 
many cases. A method requiring a vector time series 
is given by Baker, Gollub and Blackburn |ll[] and an- 
other method based on symbolic dynamics is discussed in 
Rcfs. [[T^-fTi! . The effect of noise on parameter estimation 
was studied by us Q and recently by Goodwin, Brown 
and Junge 15 1. In contrast to many of these methods our 
method in Rcf. Q works asymptotically so that an ex- 
act estimation of the parameters is in principle possible. 



The static methods based on minimization are computa- 
tionally expensive because they take a longer time to run 
due to many iterations required for convergence and they 
also require annealing to eliminate the possibility of get- 
ting trapped into a local minimum. The dynamic method 
as described in Ref. [^) requires only one time evolution 
of the system equations. The method also takes care of 
annealing in a dynamic way. 

In the first part of this paper we review our method 
for parameter estimation in brief. We then extend it to 
a case when the time series of a scalar function of phase 
space variables is given. We then go on to study the ap- 
plicability of the method to a general quadratic flow in 
three dimensions. This system has a large number of pa- 
rameters and we try to estimate some of them using our 
method. 

In the second part, we show that it is possible to ex- 
tend our method to a more realistic situation, when the 
given time series is truncated after a finite time. We 
find that a repetitive use of the finite time series can be 
made to estimate the unknown parameters of the under- 
lying system without altering the dynamic nature of the 
method. The accuracy of such an estimation increases 
with the increasing length of the given time series. We 
also see that the accuracy saturates with the number of 
times the finite time series is used. 

Lastly in the third part of this paper, we suggest an 
interesting application of parameter estimation method. 
Consider a situation where an unknown perturbation dis- 
turbs a known chaotic system. In many practical situ- 
ations when the external perturbation is unknown, an 
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ansatz function modeling the behaviour of the external 
perturbation is tried. We show that it is possible to use 
our parameter estimation method, to confirm the form of 
an ansatz function modeling the external perturbation. 

In section IIA we briefly introduce our method of pa- 
rameter estimation and discuss its important features. 
In section IIB we extend it to a general situation when 
the given time series is obtained as a scalar function of 
the phase space variables. Section IIC deals with a gen- 
eral quadratic flow in three dimensions. In section III 
we extend the method to the case of finite time series 
and present two examples. Finally in section IV we give 
the application of the method in confirming the form of 
an unknown external perturbation to a known dynamical 
system. In section V we conclude with a summary of the 
results. 



II. PARAMETER ESTIMATION 

A. The method 

Here, we briefly introduce our method for parameter 
estimation from a scalar time series. We would like to 
direct the reader to Ref. Q for a more detailed discus- 
sion. We start by considering an autonomous dynamical 
system of the form, 



x = f (x, a) 



(1) 



where x = [x\, X2, ■ ■ ■ , x n ) is an n-dimensional state vec- 
tor whose evolution is described by the function f = 
...,/„). We denote a set of m unknown scalar pa- 
rameters by a — {ot-L,a<2, . . . ,a m ). A possible appear- 
ance of any other parameters (assumed to be known) is 
not shown in Eq. (|lj). 

Without loss of generality we assume that a time se- 
ries of the variable x\ is given. The problem we consider 
is to estimate a from the given scalar time series of x\ 
assuming the functional form of f to be known. 

In analogy with the control method used earlier by 
John and Amritkar we combine synchronization 

with adaptive control to achieve our goal of estimating 
a in Eq. (|l|) as follows. We construct another system 
of variables x' having a structure identical to that of 
Eq. (|]) with a linear feedback proportional to the dif- 
ference x[ — x\ added in the evolution of the variable x\. 
Thus the system is given by, 



x\ = /i(x',a') - e(x' 1 - Xi) 

2,...,n. 



(2) 



where the function f = . . . , /„) is the same as that 
in Eq. ([!]). The initial values of parameters a' which 
correspond to the unknown parameters a in Eq. (|l]) are 
chosen randomly. The newly introduced parameter e is 
the feedback constant. It is known that if a' = a then the 
systems (Q) and (||) synchronize after an initial transient, 



provided the conditional Lyapunov exponents (CLE's) 
of the system (||) are all negative [||. The CLE's are 
obtained from the eigenvalues of the Jacobian matrix J 
whose elements are given by, 



Ji: 



dx~ - 



(3) 



Since the values a = (ai, . . . , a m ) are unknown, we 
need to set a' = (a[, . . . , a' m ) to random initial values 
and evolve them adaptively so that they converge to the 
values a. Note that a good guess for the initial values of 
a' , can be useful in many cases. 

We first consider the case when a (and its counterpart 
a') contains only a single element, i.e. the case when only 
a single parameter in Eq. ([!]) is unknown. For notational 
simplicity we now denote this single parameter by a. We 
start with a random initial value for a' and evolve it in 
a controlled fashion so that it converges to a. This is 
achieved by raising a' to the status of a variable which 
evolves as, 



-8{x' l — x\) w 



dfi 
da' 



(4) 



where S is called stiffness constant and w is some suit- 
ably chosen function of dfx/ da'. A simple choice for w 
is w = dfi/da' giving the adaptive evolution equation 
for a' as, 



a =-8{ Xl -x x ) — . 



(5) 



Eq. (Q) or Eq. (g) when coupled with Eq. (^) con- 
stitutes our method of parameter estimation. A vector 
(x', a') initially set to random values asymptotically con- 
verges to a vector (x, a) in Eq. (|l|) provided the con- 
ditional Lyapunov exponents (CLE's) for the combined 
system (Eqs. (||) and (||)) are all negative. This facilitates 
the estimation of a. 

Eq. (|^) is equivalent to a dynamic algorithm for min- 
imization of synchronization error between Eqs. (0) and 
(0) as discussed in Ref. 0. 

Note that if we assume in the above discussion that the 
unknown parameter a appears in the function fi corre- 
sponding to the variable x\ for which the time series is 
given then the calculation of the factor dfi/da' in Eq. (|5|) 
is straightforward. However this may not be necessarily 
the case. The parameter a may appear in any of the 
other system functions. If it appears in the functions for 
the variables for which the time series is not given, e.g. 
in any of the functions /2, . . . , f n in Eq. (|l|), then corre- 
spondingly the calculation of the factor dfi / da' becomes 
nontrivial. 

To make this point clear we assume that the unknown 
parameter a appears in the function fk (x) governing the 
evolution of variable Xk with k ^ 1 while the time series 
of x\ is given. In such a case Eq. (^) gets modified to, 
(See Ref. §) 
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(6) 



Further if the variable Xk itself does not appear in the 
function f\ then the complexity of the calculation still 
increases. This issue has been explained in detail with 
an example in Ref. 

Next we consider the case when the set a of un- 
known parameters contains more than one element, say 
(ai,a>2, ■ ■■)■ Now we set up an adaptive evolution for 
each of the corresponding parameters (a{, a'2, . . .). For 
the case of two unknown parameters a\ and a.2 , appear- 
ing in functions fk and /; respectively, the adaptive evo- 
lution is given by, 



dhdfk 
dx' k da 1 
s , , s dfx dfi 



(7) 



where <5i and 8% are two stiffness constants deciding the 
rates of convergence. For estimating the values of a± and 
«2 Eqs. (0) can be coupled with Eqs. (||) which provide 
the necessary synchronization of system variables if the 
associated CLE's are negative. 

In the next subsection we extend our method to a sit- 
uation when a time series of a scalar function of phase 
space variables is given. We show that it is not only possi- 
ble to build a synchronizing system but also to adaptively 
estimate an unknown parameter. 



B. Parameter estimation using time series of a scalar 
function of variables 



In our discussion of parameter estimation in earlier 
subsection, we have assumed that time series of one of 
the phase space variables is given. This may not be the 
case in many practical applications and in general the 
observed quantity can be a function of the phase space 
variables, say s(x). It is possible to construct a synchro- 
nization scheme in such a situation Jl6j . 

We consider the system given by Eq. ([!]) and assume 
that the time series s(x) which is a function of phase 
space variables is given. A synchronization scheme can 
be set up in this case by using a suitable modification of 
the feedback in Eq. (0) as follows pii] . 



9s' 



.r , = /i(x',a) - esgn ( — ) (s' - s(x)) 



(8) 



where s' = s(x') and we give a feedback proportional to 
(s' — s) in the function /i with feedback constant e. The 
function s(x) denotes the given time series. 

It can be shown that if the parameters a are as- 
sumed to be known, the above system of equations for 



x' (Eqs.j H)) converges to x, provided the CLE's are all 
negative fllqj . 

In Eqs. (H), we have assumed that s(x) has an explicit 
dependence on the variable x± so that ds'/dx[ ^ 0. If 
this is not the case, we can choose any other variable for 
the feedback on which s(x) depends explicitly. The fac- 
tor sgn() in Eq. (^) makes sure that the term provides a 
'negative feedback' for all the time so that a convergence 
is feasible. 

To estimate parameter a in such a case, we set up a 
synchronization scheme combined with an adaptive con- 
trol in analogy with Eqs. (||) and (|J). This system can 
be written as, 



/i(x',a') - e sgn 



•J-; ./Vx'.n': j = 2,. 



s(x)) 



,,n. 



" 5sgn (^[.) (S '~ S(x)) ^- 



(9) 



Eqs. ([|) can be used for estimating a when a time series 
of s(x) is given. The condition for such an estimation of 
a to be possible is that the CLE's associated with the 
system (Q) are all negative. 

To demonstrate the above procedure, we consider the 
Lorenz system given by, 

x = a(y — x) 
y = rx — y — xz 

z = xy — bz, (10) 

where the variables (x, y, z) define the state of the system 
while (er, r, b) are the three parameters. We consider the 
case when the time series of s{x,y, z) — 0.5a; 2 + l.ly is 
given as an output of the above system and the parame- 
ter a is unknown. 

To estimate the value of a, we form a system of vari- 
ables {x' , y' , z' ', er') similar to Eq. (^). The evolution equa- 
tions are 

x' = o-'(y' -x')-e sgn(z')(s' - s(x, y, z)) 



!J 



rx' — y' — x' z 



(11) 



z = x y — bz 

&' = -6 sgn(>')(s' - s(x, y, z))(y' - x'), 

where s' = 0.5x' 2 + l.ly'. 

Figures l(a)-(d) show the evolution of the differences 
x' — x, y' — y, z' — z, a' — a respectively (Eqs. ( ^p| ) and (|Tl|)) 
as a function of time t. We see that these differences all 
go to zero as t — > oo. This indicates that an unknown a 
can be estimated using Eq. (pl|). 

The CLE's are obtained using the Jacobian matrix J 
given by 

— (j — c sgn(x)x a — l.lesgn(a;) y — x 
r — z —1 —x 

y x —b 

-Ssgn(x)x(y — x) — 1.1(5 sgn(x)(y — x) 

(12) 



J = 
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We have verified that all the CLE's are less than zero 
except one trivial CLE which is zero. 

We have performed simulations and successfully esti- 
mated unknown parameters in Lorenz system with other 
forms of the function s(x,y, z). The function s(x,y,z) 
should however be such that all the associated condi- 
tional Lyapunov exponents should be negative. 

C. A general quadratic flow in 3-D 

Now we consider a quadratic flow in 3-D given by, 

x = ao + a±x + a 2 y + a^z + a^x 2 + a^y 2 

+ a,QZ 2 + a 7 xy + a$yz + agxz 
y = b + bix + b 2 y + b 3 z + 6 4 .x 2 + b 5 y 2 

+ bftZ 2 + b'jxy + b%yz + bgxz 
z = c + c\x + c 2 y + C3 z + C4X 2 + C5?/ 2 

+ cqz 2 + c-jxy + c$yz + cgxz, (13) 

where (ao, . . . , ag, bo, . . . , bg, Co, . . . , eg) form a thirty di- 
mensional parameter space and [x, y, z) are the three 
variables. We have performed simulations in which we 
have assumed more than one of the thirty parameters 
of the system ( |l3|) to be unknown and tried to estimate 
them when a time series of one of the variables is given. 

To elaborate, we assume some of the thirty parameters 
to be unknown while the remaining to be known. Some of 
the known or unknown parameters may be zero thereby 
making the corresponding term absent from the system. 
To illustrate the procedure we consider a case when three 
parameters (01, a 2 , ay) are unknown and a time series of 
x is given, we set up a system of equations similar to 
Eq. (||) with the adaptive control loops similar to Eq. (Q) 
for the three parameters (a'i,a' 2 ,a 7 ) as, 

dj = — Si(x' — x)x' 
a' 2 = -S 2 (x' - x)y' 

a' 7 = -S 3 (x' - x)x'y'. (14) 

Eqs. ( |l4| ) when coupled to the system of variables 
(x',y',z') with an identical structure of evolution as 
Eq. ([l3]) with a feedback term in the evolution of x', can 
provide the necessary estimation of parameters when the 
CLE's associated with the reconstructed system are all 
negative. 

In Fig.2(a)-(c) we plot the time evolution of the differ- 
ences a[ — Cbi, a'2 — ct2) a' 7 — aj as a function of time. The 
correct value of a 7 was zero while the other two were non- 
zero. All the differences go to zero indicating the feasi- 
bility of simultaneous estimation of the three parameters 
(01, fl2, a 7 ) even when the actual value of one of them is 
zero. This shows that the method does not falsely detect 
a term which is absent in the system. 

We have found cases when our method can be used 
successfully for the system (|l3]) to simultaneously esti- 
mate as many as five parameters. (One such case is the 



set of parameters ai, a 2 , a 7l 63, c\, while the time series of 
x is given.) 

Further we have also found that when any two of the 
thirty parameters in the system ( |l3| ) are unknown, we 
can apply our method to simultaneously estimate them 
asymptotically to any desired accuracy when the time se- 
ries of a suitably chosen variables is given. Our results 
suggest that the information about all the thirty param- 
eters should in principle be contained in the time series 
of a single variable of the system, though at present we 
do not have any systematic approach to the simultaneous 
estimation of all of them. 



III. PARAMETER ESTIMATION USING A 
FINITE TIME SERIES 

A. Algorithm for repetitive use 

In this section we discuss an algorithm for repetitive 
use of our method to impove the accuracy of parameter 
estimation when the given time series is of finite duration. 

Before going on to describe the algorithm it should be 
mentioned here that even if a finite time series is used 
repeatedly, we do not expect an exact estimation of the 
unknown parameter. A finite chaotic trajectory sets a 
limit on the accuracy to which the unknown parameter 
can be estimated. This can be seen as follows : 

We consider symbolic dynamics on the attractor which 
provides a generating partion of the attractor. It is well 
known that as the system evolves in time, a finer and finer 
coarse graining is required to specify a particular trajec- 
tory or alternatively, the trajectory gives us a finer coarse 
grained information about the attractor. The number of 
coarse grained partitions as a function of time goes as, 

n p ~ exp{ht}, (15) 

where h is the Kolmogorov entropy fl7|| . 

If £ d is the volume of a hypercube in a d dimensional 
phase space and if the size of the attractor is normailized 
to unity, the number of hypercubes in a generating par- 
tion may be approximated as, 

n p ~ (16) 

Equations ([l5]) and ([l6]) indicate that the length scale of 
a hypercube in a generating partition goes as, 

£ ~ exp{-- d t}. (17) 

It can be seen from Eq. ( |l7| ) that as long as t is finite, 
the volume of the hypercube in a coarse graining of the 
attractor will not reduce to zero. Thus a finite trajectory 
sets a limit on the accuracy to which any information 
can be extracted from it. This can be further related 
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to Lyapunov exponents using the famous Kaplan- Yorke 
conjecture |Q as, 

£ ~ exp{-^^t}, (18) 
a 

where A is the characteristic Lyapunov exponent of the 
system. For a chaotic system with a single positive Lya- 
punov exponent denoted by A + , eqn. ( |l8|} reduces to, 

£ ~ exp{-^-t}. (19) 

Now we will discuss the algorithm for repetitive use of 
a hnite time series to estimate an unknown parameter. 
Similar to the case considered in section IIA, we assume 
that the parameters a = (ax, ct2, ■ ■ ■ , ot m ) in Eq.(|l|) are 
unknown while the time series of x\ is given. We further 
assume that the time series is truncated after a finite time 
T. 

For the time interval < t < T, we can use the proce- 
dure identical to that described earlier (Eqs. (|J) and (|^)) 
to evolve variables (x', a') with random initial conditions. 
The given finite time series is fed in system ^ as in the 
earlier case. In this way we can get an approximate value 
of a which we denote as, a 1 = a'(T). 

Now at time t = T we set the variables x' to exactly 
the same (randomly chosen earlier) initial values while 
a' = a 1 and feed the same finite time series {x(t)\0 < 
t < T} again into the system (|J) through the feedback 
terms in Eqs. (||) and (||), i.e. we set x(t + T) = x(t). 
We now evolve the variables (x', a') for the time interval 
T < t < 2T to obtain a new estimated value of a which 
is a 2 = a'(2T), 

We repeat the procedure to get successive estimates for 
the value of a denoted by a 1 , a 2 , a 3 , ... , a N , ... at times 
t = T, 2T, 3T, . . . , NT, . . . respectively. Thus, starting 
from an initial guess for the value of a we obtain a se- 
quence of estimates a , a 1 , ... , a" after N usages of the 
given finite time series. For large enough N we get a 
better and better estimate of a, although eventually the 
accuracy of such an estimate saturates as N is increased 
further. 

The conditions for the method of parameter estimation 
using a finite time series to work successfully are : 

1. The conditional Lyapunov exponents associated 
with the reconstructed system should be all neg- 
ative. 

2. The time T after which the given time series is trun- 
cated should satisfy T > r where r denotes the 
transient time required for synchronization of the 
systems (Q) and (||) with the parameter evolution 
given by Eq. (||). 

In the next subsection, we discuss two examples of pa- 
rameter estimation from a finite time series, viz. Lorenz 
system and an electrical circuit of a phase converter. 



B. Examples 

1. Lorenz system 

As our first example we choose the Lorenz system 
given by Eq. |[c]) where we assume that the time se- 
ries {x(t)\0 < t < T} is given and the value of a is to be 
estimated. We set up the following system of equations 
(see Eqs. (|) and (|)). 

x = a(y' — x) — e(x — x) 

■i i i ii 

y = rx — y — x z 

z' = x'y' - bz' 

&' = -8(x' - x)(y' - x') (20) 

where we feed the given time series in the evolution of x 
for the interval < t < T to obtain the first estimate a 1 . 

As described in the earlier subsection, we then go on 
repetitively feeding the same finite time series x(t) in 
Eq.(p0[) to obtain successive estimates for the value of 
a. Starting from a random initial value we denote this 
sequence of estimates by a , a 1 , ... , cr N where N denotes 
the number of times we use the given time series. 

In Fig. 3 we plot the evolution of the difference a' — a as 
a function of time t during the time interval < t < 3T 
where we use the time series x(t) thrice. We see that the 
difference decreases as we increase the number of times 
the finite time series is used. We also observe that shortly 
after each resetting of the initial vector (x\y', z') which 
is done at times T, 2T, the synchronization weakens and 
fluctuations are present. This is due to the random reset- 
ting of the y and z components which gives a transient 
before the synchronization is recovered. An appropriate 
feedback constant e may be chosen to lessen this transient 
in every usage of the time series. 

In Fig. 4 we plot the successive differences <j n — a as a 
function of AT, the number of times we use the given finite 
time series. We see that the difference a N — a goes on de- 
creasing with increasing N. However, as N is increased 
further, it saturates to a constant finite value depending 
on the length of the time series used for the calculations. 
This is consistent with our expectations that finite time 
series can contain only finite information about the sys- 
tem as discussed in the previous subsection, e.g. using 
A + ~ 0.9, the finest length scale that can be obtained 
using a finite time series with T — 30 is estimated to be 
0.05 (Eq. (|l~9l)) which means an accuracy of about 10~ 3 . 
This is also the order of magnitude of the accuracy of 
parameter estimation. 

The three curves in Fig. 4 correspond to three different 
values of T = T\ < Ti < T3. We see that an increasing 
T gives better estimate of the parameter. This is natural 
since a very long time series corresponding to T — > 00 
is expected to give an exact estimation of the unknown 
parameter. 
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We have similarly implemented our method to estimate 
other parameters of Lorenz system using finite time series 
of either x or y. The method fails to estimate any of the 
parameters when the time series of z is given. The rea- 
son for this is that one of the associated conditional Lya- 
punov exponents is critically zero and the convergences 
are slow. 



2. A phase converter circuit 

As our next example, we consider the set of equations 
describing an electrical circuit for a phase converter Jlfitl 
system in a dimensionless form given by, 

±1 = x 2 

x 2 = —KX2 — — (x 1 + 3x 3 ) 

X 3 = Xi 

x A = -kx A - ^ (x\ + 3x\) + Boost (21) 

where k and B are the two parameters. Here we consider 
the time series {x2(i)|0 < t < T} to be given. Notice that 
the system has a simple time dependent term making 
it a non-autonomous system. Such a system is equivalent 
to an autonomous system in higher dimensions. We have 
successfully estimated any one of the parameters k or B 
(or both) using finite time series of X2 (t) . 

Figure 5(a) shows a schematic diagram of the circuit 
for the phase converter. The system in known to exhibit 
a chaotic behaviour due to period doubling bifurcations, 
codimension two bifurcations etc. Figure 5 (b) shows a 
chaotic attractor in the x\ — X2 plane of the phase space. 

Figure 6 shows the plot of the successive differences 
k N — k as a function of N, the number of times we use 
the given time series for two different values of the trunca- 
tion time T. As expected the accuracy of the estimation 
increases with increasing T while showing a saturation 
with increasing number of repeated usages. 

Thus, we have shown how the method of parameter 
estimation can be used when a finite time series is given. 
The method works when the CLE's associated are all 
negative and the time series given is of longer duration 
than the transient time required for synchronization. 

IV. FORM OF A MODEL PERTURBATION 

Here we describe an interesting application of our 
method to test a function modeling an unknown exter- 
nal source of perturbation to a known chaotic system. In 
many practical situations when an external source of dis- 
turbance in not known, a trial function is used to model 
the perturbation. 

We imagine a situation when it is required to verify a 
proposed trial model form for the perturbation. We de- 
note the actual perturbation by a function F(x, jj) and 



the trial function by G(x', //) where /x and /x' are param- 
eters. In the following, we demonstrate the use of our 
method of parameter estimation to confirm the form of 
the trial function. Note that here we do not deal with 
the issue of obtaining the form of the model function. 

Now if the proposed trial function G models the ex- 
ternal perturbation F correctly, then a scheme based on 
synchronization combined with adaptive control should 
produce synchronization of variables and make the pa- 
rameters // converge (to /i). Thus a successful synchro- 
nization should then indicate a correctly chosen model 
function. In this manner we can use the method to dis- 
tinguish between a correct model and a wrong model for 
an external perturbation. We elaborate on this applica- 
tion further using the example of Lorenz system. 

Consider the Lorenz system perturbed by a sinusoidal 
term F — Asin(wx), 

x = a(y — x) + Asm{u>x) 

y = rx — y — xz 

z = xy- bz, (22) 

where we assume the unperturbed Lorenz system to be 
known. The function F — Asm(ujx) is the external per- 
turbation. We assume that the time series of x is given 
as an output of the system (|22]). 

To set up the required scheme we construct a system 
of variables {x 1 , y' , z') and their evolution as, 

x = a(y' - x') + G(x', y', z', //) - e(x' - x) 
y = rx' — y' — x z 
z = x y — bz', 

(,' = -5{x' - x) — . (23) 

where G(x' ,y' z') is the trial perturbation function. 

We feed the time series x(t) obtained from system (|2^ ) 
into the model system (|23|). Now if G models the be- 
haviour of F correctly then the two systems should ex- 
hibit synchronization while the parameters should show 
convergence to the correct values. In our simulations we 
have tried several different forms for the trial function G. 

Figures 7(a)-(c) show the time evolution of x' — x, fii 
and [ii respectively while the feedback is given into x and 
the trial function is G — \i\x 2 +/i2- It can be clearly seen 
that there is no synchronization of variables. The trial 
function G = ^ix 2 + [i% thus fails to produce synchro- 
nization and hence can be discarded as a plausible model 
for F. We also note that the parameters fi[ and ^' 2 do 
not show convergence. 

In Figs. 8 and 9, we plot similar graphs for two more 
choices of the trial function. In Fig. 8(a)-(c) we use 
G = fiix — fi2X 3 and plot x' — x, [i\ and /i2 respectively. 
We choose this form of G since it represents the two 
leading terms in the series expansion of the the func- 
tion F — Asm{ijjx). We can see from Fig. 8 that such an 
approximation fails to produce synchronization and also 
the convergence of parameters. 



G 



As a third choice we use G — fix sm((X2%) in Eq. ( p3] ) 
and plot the time evolution of x' — x, /ii and [i% in 
Fig. 9 (a)- (c) respectively. The difference x' — x goes to 
zero as time increases showing synchronization. The pa- 
rameters \i\ and \ii converge to the correct values A and 
u) respectively. The variables y' and z' also synchronize 
with y and z respectively. This confirms that this trial 
function correctly models the function F. 

Now as a last consideration, we use the form G = 
Hi sm(/i2x) again but unlike in Eq. ( ^3|) we perturb a 
wrong variable in the model system, i.e. we choose to 
add the trial perturbation in the evolution of say y' . The 
feedback is given in x. The evolution equations are 

x' = a(y' — x 1 ) — e(x' — x) 

y' = rx'-y'-x'z' + G(x',y',z',fj,') 

z = x'y' — bz' , 

c>C 

= -5{x' -x) — . (24) 

In Fig. 10 (a)-(c) we plot the time evolution of x' — x, [i\ 
and /x 2 respectively. We see that even if G correctly mod- 
els F, synchronization does not take place. This shows 
that along with the form of F we can also confirm a guess 
about the perturbed variable. 

Thus, the results presented in this section suggest that 
the method which we use for estimating parameters can 
be used to distinguish between a correct trial function 
and the wrong trial functions for an unknown external 
perturbation to a known system |2(i|j . 



after a finite time. We have shown that repetitive use of 
a finite time series can be made to estimate an unknown 
parameter of the system. The accuracy of the parameter 
estimation saturates as the given finite time series is used 
more and more number of times. The accuracy increases 
with the increasing length of the given time series. 

In the end we have demonstrated an important applica- 
tion of our method in confirming the correctness of a trial 
model function for an unknown external perturbation to 
a known system. We see that a perfect synchronization 
between a perturbed system and its dynamical copy us- 
ing a model for the perturbation is possible only when 
the form of the trial function is correctly guessed. These 
results indicate that our method can be used as a test for 
the trial model for an unknown external perturbation to 
a known system. Another possible application (not dis- 
cussed in the paper) is as follows. Our method may be 
employed to experimentally measure the unknown value 
of a component added to a known circuit. In such a situ- 
ation the equations governing the circuit are known, and 
can be used to estimate the unknown component value 
accurately. This is feasible due to the asymptotic con- 
vergences in our method. 
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V. SUMMARY AND CONCLUSIONS 



We have described dynamic method of parameter esti- 
mation from a given chaotic time series of a phase space 
variable of a dynamical system 0. Further, We have 
generalized the method for the case when the quantity 
for which the time series is given is a scalar function of 
the phase space variables. We have shown that it is not 
only possible to synchronize two systems using the time 
series of the scalar function but also to asymptotically 
estimate unknown parameters adaptively to any desired 
accuracy. This is done by providing a linear feedback in 
the evolution of one of the variables on which the scalar 
function explicitly depends. The method works success- 
fully provided the function for which the time series is 
given is such that the associated conditional Lyapunov 
exponents are all negative. 

We have also applied our method to a system with a 
large number of parameters, i.e. a general quadratic flow 
in 3-D. We have observed that a simultaneous estimation 
of a few parameters is possible provided the condition of 
convergence as stated in Ref. || is satisfied i.e. all the 
CLE's are negative. 

As a next consideration, we have extended our method 
to a realistic situation when the given series is truncated 
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Fig.l (AM &c REA) 



FIG. 1. The plots (a)-(d) show the evolution of the differences x' — x,y' — y,z — z,a' — a as a function of time for the 
Lorenz system (Eqs. ( |Io| ) and (p"l|)) respectively for the case when a time series for s(x,y,z) — 0.5a; 2 + l.ly is given and a is 
unknown. The differences go to zero asymptotically indicating that it is possible to use our method to estimate an unknown 
parameter when time series for s(x) is given. 
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(AM & REA) 



FIG. 2. The plots (a)-(c) show the evolution of the differences a[ — ai, 02 — 0,2, a'r — 0,7 for a general quadratic flow in 3-D 
(Eq. (|l|) and ([ij)), plotted as a function of time when the time series of x is given. We see that all the differences approach 
zero indicating the feasibility of simultaneous estimation of more than one parameter. The correct value of a-j was zero, showing 
that a term absent in the flow equations is not falsely detected by our method. 
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Fig. 3 (AM SuBEA) 

FIG. 3. The plot shows the evolution of the difference a' — a as a function of time in the Lorenz system (Eq. (|20|)) with 
unknown parameter a when the given time series of x is truncated after the time T = 20. We have used The finite time series 
thrice and plotted the curve for the interval < t < 3T. We see that the successive values of the difference at t = 0, T, 2T, 3T 
decrease. This indicates that a repetitive use of the finite time series can improve the accuracy of parameter estimation. 
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Fig.4 (AM & REA) 

FIG. 4. The graph shows the successive differences a N — a plotted as a function of N, the number of times a finite time 
series {x(t)\o <t< T} is used to estimate an unknown a in a Lorenz system (Eq. (^o|)). We see that after an initial transient, 
the difference decreases showing better accuracy of the estimation. We also see that as N increases further the accuracy of 
estimation saturates and it is not possible to improve upon the estimation beyond this. The three curves correspond to three 
different values of T where Ti < T2 < T3. It can be seen that a larger T leads to a better estimation as expected. 
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Fig.5(bi (AM & REA) 

uitlEq. (BID) which shows a chaotic behaviour. Fie. 5(b) 



FIG. 5. A schematic diagram (Fig. 5(a)) of a phase converter circuit^Eq. ([2l|)) which shows a chaotic behaviour. Fig. 5(b) 
shows a chaotic attractor for the parameter values k = 0.1, B = 3.0. 



14 



1.5 




20 40 60 80 100 

N 



Fig.6 (AM & REA) 



FIG. 6. The graph shows the successive differences k — k plotted as a function of N, the number of times a finite time 
series {x2(t)\o < t < T} is used to estimate an unknown A: in a phase converter circuit system (Eq. (|2l|)). We see that after 
an initial transient, the difference decreases showing better accuracy of the estimation. We also see that as N increases further 
the accuracy of estimation saturates and it is not possible to improve upon the estimation beyond this using our method. 
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Fig. 7 (AM &c REA) 



FIG. 7. The plots (a)-(c) show the time evolution of x' — x, fii and \i2 respectively for the Lorenz system with the feedback 
given in equation for x and with the trial perturbation function G = fiix 2 + fj,2 while the correct perturbation is F = Asin(u>x) 
(Eq. (|23|). We see that the guess function G — fii(Eq. (p3|))a: 2 + ^2 fails to produce synchronization and hence can be discarded 
as a plausible model for F. It can also be seen that there is no convergence of the parameters taking place. 
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Fig. 8 (AM & REA) 

FIG. 8. The plots(a)-(c) show the time evolution of x' — x,(Mi and fi2 respectively for the Lorenz system with the feedback 
given in equation for x and with the trial perturbation function G = (i\X + H2X 3 while the correct perturbation is F = Asin(u>x) 
(Eq. ([23])). It can be clearly seen that even when G = fiix — fi2X 3 matches in form with F upto two leading terms in the 
expansion of F, it fails to produce synchronization and hence can be discarded as a plausible model for F. Also there is no 
convergence of the parameters taking place. 
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Fig. 9 (AM &c REA 

FIG. 9. The plots (a)-(c) show the time evolution of x' — x, fii and /12 respectively for the Lorenz system with the feedback 
given in equation for x and with the trial perturbation function G — fii sin(fi2x) while the correct perturbation is F — Asin(uix) 
(Eq. (^3|)). It can be clearly seen that the difference x' — x converges to zero asymptotically indicating an exact synchronization 
between the variables. Thus by using our method the guess for the model perturbation function can be easily justified. 
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Fig. 10 (AM & REA) 



FIG. 10. The plots (a)-(c) show the time evolution of the difference x' — x and the parameters fii and fi2 respectively for 
the Lorenz system with the feedback given in the equation for x and with the trial perturbation function G = /Lti sin(/^2^) in 
the equation for y while the correct perturbation is F = Asm(ux) in the equation for x (Eq. (p4|)). Thus, unlike the case 
plotted in Fig. 9 the trial function used here, perturbs the wrong variable. It can be clearly seen that the trial function G does 
not produce synchronization between variables. The parameters also do not converge. Thus as expected, the guess function 
G — /ii(sin fi2x) when added to a wrong variable, cannot model the perturbation. 
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